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1. Introduction 

The Hhota's direct method (hereafter the direct method) provides us with a very 
powerful tool in the integrable systems Q. Nakamura applied the direct method to 
the Ernst equation and obtained Tomimatsu-Sato solution(TS solution) in bilinear 
forms 1^. However, his bilinear form does not take completely the same form as the 
conventional bilinear forms in the following senses. It can not be expressed only 
by the Hirota's derivatives but involves ordinary derivatives. Also the coefficients 
of the Hirota's derivatives are not constant but functions of independent variables. 
Therefore, it was not trivial that the direct method works well in this system. In 
the previous paper, we proved that the direct method does work in this system^. 
However, our proof was complete in the restricted one dimensional case, Weyl solution, 
and was incomplete in full two dimensional case, TS solution. Naive Pfafhan identity, 
which was valid for one dimensional case can not be applicable to double Wronskian 
in two dimensional case. We consider the origin of this trouble lies in the peculiarities 
of the bilinear form mentioned above. By adopting the multilinear form|^, |^, we can 
rewrite the above bilinear forms so as to involve only multilinear operators. Thus we 
are forced to go beyond the bilinear form. However, so far, any trilinear equations have 
not been shown to be integrable explicitly. In this paper, we prove the integrability 
of the Bogoyavlenskii-Schiff(BS) equation^, |7[ Furthermore, being guided by the 
strong symmetry we search an integrable system in (3 + 1) dimensions. 

This paper is organized as follows. In Sec. 2, we construct the exact N soliton 
solution of the BS equation in N x N Wronskian representation. In Sec. 3, a 
constructive proof of the N soliton solution is given from the Miura transformation 
and the Hirota condition. In Sec. 4, we propose a new equation in (3 + 1) dimensions 
by the strong symmetry and give the travelling solitary wave solution to this system. 
Sec. 5 is devoted to discussions. 

2. Exact N Soliton Solution of the BS Equation in N x N Wronskian 
Representation 

We review the treatment to find the exact solutions of the KdV equation in the 
direct method for later use. The KdV equation is written as 



Ut + <i>{u)ux = 0, (1) 

where ^{u){= + + 2uxd~^) is the strong symmetry^. Potential form of this 
equation is 



(m = <f)x)- 



(2) 



By the dependent variable transformation 




(3) 



T 




(4) 



where the Hirota's derivative T) operating on / ■ 5 is defined by 



V:J{z) ■ g{z) ^ {d.. 



(5) 



We have, in general, an exact solution r^y which can be expressed as 
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N 

Tat = 1 + ^ ^ exp(Ai^ H h Ai„), (6) 

n=l IV C„ 



{Po - Pk 



l2 



where AfC^ indicates summation over all possible combinations of n elements taken 
from N, and symbols Cj always denote arbitrary constants. Equation (6) together 
with u = 2 (log r) 2:3; gives N soliton solution of the KdV equation [Q. 

Then we proceed to the study of the BS equation which can be described not by 
the bilinear form but by the trilinear form. The BS equation is given by 

ut + <^{u)u^ = 0, (10) 

Here $(w) has the same form as that in equation (|l|) with argument x. Using the 
potential u = (j)x, this equation reads 

4>xt + 4>xxxz + '^i'xi'xz + '^4>xx4>z = 0. (11) 

This equation has been constructed by Bogoyavlenskii and Schiff in the different ways. 
Namely Bogoyavlenskii used the modified Lax formalism]^, |^, whereas Schiff obtained 
the same equation by the reduction of the self dual Yang-Mills equation^ . In reference 
[Q, ||, ^ , it was shown that equation ( pi] ) is transformed into the trilinear form 

%{r^T: + 8r,2r;r, + 9%%)t • t • t - o, (12) 

through the dependent variable transformation(|^). The operators T, T* are defined 

byii 

r,"/(z) -5(2) • h{z) = +jd,, +fd,,rf{z,)g{z2)h{z3)U,=,,=,s=,, (13) 

where j is the cubic root of unity, j — exp(2i7r/3). T* is the complex conjugate 
operator of % obtained by replacing {dz^ + jdz2 + j^dz^) by (9^^ + j^dz2 + jdz^). To 
find the N soliton solutions, we repeat the same procedure as in the case of the KdV 
equation. We find that tn is expressed as 

N 

TAf = 1 + ^ ^ 77ii...j„ exp(Aii H h Ai„), (14) 



where 



Xj = PjX + qjZ + Tjt + Cj, Tj = -PjQj- (15) 



The proof is referred to Sec. 3. In the case of = 2, the above 2-soliton solution is 
same as that obtained by Schiff S . 
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We rewrite tn in the form of N x N Wronskian, 



tn — det 



/ /i 



N 



"x In , 



where 



!i = exp 



-(pjX + gjZ + r^i + Cj) 



exp 



(16) 



(17) 



The degree of variables in typical soliton equations are fixed. For example, the KdV 
equation (0) demands that 



3 [a.] - \dt\ 

where is the degree of d^- So we may set = 1, 
\d^\ = 1, \dt\ = 3. 



(18) 



(19) 



We can use the Wronskian technique for the Wronskian solutions of the KdV 
equation||l^. However, it is not the case in the BS equation. Since equation ( p^ 
only demands 



equation (^0|) allows an indefinite factor, say a, like 
[9,] = 1, [9,] = a, \dt\ = 2 + «. 



(20) 



(21) 



In this case, we cannot use the Wronskian technique by the presence of an indefinite 
factor a. This may enforce us to extended the Pfaffian identities. We checked that 
( p^ ) are solutions to equation (|l^) for an arbitrary a by the computer program 
Mathematica to = 8. 

Figure |l| shows an example of the propagation of one soliton(M). The potential((/>) 
corresponding to Figure |l] with two floors is shown in Figure ^j. In Figure ^, typical 
patterns of two solitons(pi ^ pi) and the potential with four floors are depicted. 
In the soliton collision, however, appears a new feature. For the special momentum 
combination (pi = P2 7^ 0) two solitons shrink to V form(Figure |^): we may call this 
pattern V soliton. 

V soliton is a peculiar feature of the BS equation. So let us discuss about it in 
more detail. In the KP equation 



/ X 

the resonance condition, 

^(kg) =w(k]_)±w(k2). 



(22) 



(23) 



and 



kg ^ ki ±k2. 



(24) 



with = {Pj,qj) gives [0, 

{pi ±p2f-A{pi ±P2){UJI ±LJ2)+3iqi±q2f = ±3piP2 ( {pi ±P2? ~ {h -hf\=^. (25) 



where Ij = qj /pj- Here r2 = 1 + eJ^^ + e^^ + rii2e^'^^^'^ with \j — pjX + qj-y + uijt + Cj 
{p'j — 'ipjiLJj + — 0) and the phase shift ryi2 is 

^ jPi ~'P2)^ -Mpi -P2)ioJi -a'2) +3(gi -q2f ^ jpi -P2f - (^1 -I2? .2fi^ 

'^'^ bl+P2)4-4(pi+P2)K+C^2)+ 3(91 +92)2 (pi+p2)2-(;i-;2)2- ^ ^ 

So the resonance condition corresponds to 7712 = or 00. In the BS equation the 
resonance condition ( p^ ) gives 

h _ P2± 2pi 

Zi~^Pi±2p2' ^^^^ 



and ?7i2 is 



.12 = (^)^ (28) 

Pi +P2 



Thus the resonance condition corresponds to neither 7712 = nor 7712 = 00. To 7712 = 
corresponds V sohton. The sohton properties of V sohton are seen from the cohision 
process of two V soHtons. Two V sohtons suffer phase shifting but conserve their 
sohtary forms after colhsion. 

3. Analytical Proof of N Soliton Solutions to the BS Equation 

We give the analytical proof that equation(^4|) is the solution to the BS equation 
(|rT|). Firstly we introduce the modified Bogoyavlenskii-Schiff(mBS) equation which 
is deduced from the Miura transformationQ. This transformation connects the BS 
solution with the mBS solution. The mBS equation is described by the coupled 
bilinear forms and tractable in the conventional Direct method. Nextly we prove 
the integrability of the mBS equation. This complete the proof of the BS solution. 

Now we proceed to the concrete explanations. We perform the Miura 
transformation in the dependent variable of the BS cquation(pT|) . 



..2 



= V 



{a = ±1). (29) 



Then we obtain the mBS equation, 

vt - Av^v, - 2v,d-^{v^), + f,,, = 0. (30) 

Equation(pO|) is reduced to the modified KdV equation in the case of a; = z. 
Introducing the new dependent variable ■0 by w = i/'a; (^), equation(|30|) is reduced to 
the potential mBS equation 

V't - 27/^.5-1 (V-^), + Vx.z = 0. (31) 

In order to eliminate the operator d^^ we describe this equation in terms of the coupled 
system, 
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p,, + = 0, (32) 

-0t + 2lpxp^z + IpzPxx + iplipz + i>xxz = 0. (33) 

Eliminating p, it is easily checked that equation(32) and equation(33) are equivalent 
to equation (|3^). Here we perform the transformation of the dependent variables, 

V-^logf^Y (34) 



p = log(FG), (35) 
then equations(32), (33) are reduced to the bilinear form, 

VlF • G = 0, (36) 
(Vt + VlV,)F -G ^0. (37) 

N soliton solutions of equations(36), (37), which we denote Fn, Gjy are speculated 
from the conventional Hirota's Direct Method, 

N 



Fn = l + J2Yl ^'^i-*" exp(A,, + ■ • • + A,J, (38) 

n=l jvC„ 
N 

Gn = 1 + J2Y1 (-l)"^*!- '-. exp(A,, + • ■ • + A,J, (39) 



n=l jvC„ 



where F/v is the same N soliton solution of the BS equation(|l4[). The proof is due 
to the Hirota condition We can rewrite the bilinear mBS equations(36), (37) as 
follows. 



^IfN • = 0, (40) 
{Vt+VlV,)fN-rN^O. (41) 

Here 



N , N N ^ 

In = ^ cxp f ^ pj (Aj + i-) + ^ PjPkAjk j , 

M=0,1 = l l<j<k ' 

fN= XI ^^'(•^J " * 9) + H ^j^kAjk I , (43) 



(42) 



£=0,1 

»p(.,.,.,„, = |a^. (44, 

denote the summation of pj — 0, 1, vj = 0, l(j = 1, 2, • • • , N). Substitution 

of the expression for and into equations (40), (41) reveals that the coefficients 
of exp(^" Xj + X]jln+i ^-^i) vanished for the respective n and m, 
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N N / , N x2\/JV. N V 

Y\\Y^^^l^ "ip^ ) r'^'P ( H y (''j - ) + H + ''j ''^ ^^j''^ ) 

^=04^=04 \ ^j=l ''^ / ^j=l l<j<fe ^ 

X cond(^, = 0, (45) 

N N I ,N \^ / ^ \ / ^ \\ 

H 51 ( Y^^^i ~ "i^P:! ) ( H^^J " ''j)* ) " ( " ''jH^^J ) 

;^=o,i £=0,1 \ ^j=i ^ ^i=i '^^ ^j=i ^/ 

/ N . N s 

X expf ^ y(^j - i/j) + (MjA^fc + i^ji^fc)^jfc j coiid(^, j/)n„t = 0, (46) 



where 



1 , for j = 1 , • • • , n : /Xj + ia,- = 1 , < n < TV 



, / N I 1 , for 7 = n + 1 , • • • , m : lii = ^ 1, n < m < N , 

— ' IX, tor J = m + 1 Jy : fij = i^j ^ ^ ' 

0, otherwise 

Here we have used the notations by J. Ablowitz and H. Segur |]l3|. As is easily seen, 
the first case of cond(/i, i£)nm gives the non trivial contribution. Equations(45) and 
(46) are reduced to the following equations (48) and (49), respectively for a given n. 



E ( E "^P^ ) ( y E '^0 n (''iP^ - ^^p^)^ = 0' (48) 

CT = ±1 V S = l ^ / ^ j = l ^ j<k 

n / / n \ 2 / n \ / n \ \ / . n 

E ( E '^ip^ ) (H'^oio) ~ [Y. ^iph^ ) °^p ( y E ^1 

n 

X n^'^J^J ~ cTfcPfc)^ = 0, (49) 
j<k 

where 

(7j=Aij-i^j. (50) 

Equation(48) is easily verified for n = 1,2. Let as denote the left-hand side of 
equation(48) as A(n) . Then A(n) has the following properties: (i) A(n) is a symmetric 
homogeneous polynomial of pj, (ii) if pi = then A(n) = 0, (iii) if pi = p2 then 

n 

A{n)^Apll[{pl-pirA{n-2). (51) 

fe=3 

Now we assume that equation(48) holds for n — 2. Then, using the properties (i), (ii) 
and (iii), we find that A{n) can be factored by a symmetric homogeneous polynomial 

f[ps n ipi-pir^ (52) 

J=l l<J<k 
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of degree n^. On the other hand, equation(48) shows the degree of A(n) to be n'^—n+2. 
Hence, A(n) must be zero for n. 

Next we discuss equation(49). we can rewrite equation(49) as 

n 

5^(7,A,(n) = 0, (53) 

i=i 

where equation (|5^) is a symmetric homogeneous polynomial of {pj,qj). 



I[(pI+p'',))(J2'^^pA e^pfyE'^O n {'^m-wf]. (54) 



etc. The first term of the right-hand side of equation(54) must be zero because this 
term contains only the odd powers of each Uj{j = 2, ■ ■ ■ ,n), the second term equal to 
zero from equation(48). Hence, equation(49) holds. 

Therefore equation(^4|) is the soliton solution of the BS equation from the Miura 
transformation(p9[) . This completes the proof. 



4. A New Equation in (3 + 1) Dimensions and Its Travelling Solitary 
Wave Solutions 

We have studied how the KdV equation in (1 + 1) dimensions is extended to the KP 
equation and the BS equation in (2 + 1) dimensions. Namely, we have two different 
ways to the integrable systems in one higher dimensions. So further analogy leads us 
to the new systems in two higher dimensions, (3 + 1) dimensions, 



'Auf + $(w)w2 ) + 3uyy = 0. (55) 
These extension schemes are schematically written in the following form: 



KdV equation (|l|) =^ BS equation (|T 
KP equation (p3) => Equation (pq) 



Equation ( |55| ) was expected to be integrable. However, the potential form of equation 
(11), 



Et + (pxxxz + '^4>x4'xz + 2(j)xx4'z + 54>yy =0, (u = (px) 



(56) 
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has a movable logarithmic branch point in the sense of WTC method|14[ . Furthermore, 
we can not construct iV(> 2) soliton solution of trilinear form of equation ( p5| ) 

(r,4r; + st^t^t, - mr^% + Ti%r^)T ■t-t = o, (57) 

by the direct method. We require the existance of 2 soliton solution. If 2 soliton 
solution, 



T2 = 1 + exp(Ai) + exp(A2)77i2 exp(Ai + A2), 
Xj = pjX + qjy + VjZ + Sjt + Cj . 



(58) 
(59) 



then 



aplplipi +P2Y - {qiP2 - q2Pi) 
ri = api, r2 = ap2, (61) 



"^'l^'2Vi^l ^ F2) ~ \H^F2 — H2F1) ,„„s 
= ..^2^2(„, ^„„^2_/.,„„_„.^ ,^2 ' (60) 



where a is arbitrary constant, thus equation(|57|) is reduced to (2 + 1) dimensional 
equation. These suggest that equation ( |55| ) is not integrable. However, equation (p5|) 
has explicit travelling solitary wave solution by tanh- function method(TFM) jl5|. The 
ansatz is expressible as a polynomial in terms of a tanh function, so that it has the 
form 

M 

u{x,y,z,t)^U(ri)^^a,T\ T = tanh(fc?7), (62) 

1=0 

where rj = x + ly + mz — ct + constant. Substitution of equation ( |6^ ) into equation 
(pa) yields an ordinary differential equation for U{i]) 



(4c + ■iV)U + 3mC/^ + m-— = 5, (63) 



where h is an integrable constant. 

We balance the highest power of T in the second term in equation (|6^) with the 
highest power of T in the final term in equation ( |63| ) to obtain 2M = Af + 2, so 
that M ~ 2. In order to solve equation (|6^) we use the automated tanh- function 
mcthod(ATFM) JTBt , where one inputs the commands in Mathematical and obtain the 
outputs in the following ways: 

In[l] := « atfm' 

In[2]:= neweq = (4 c + 3 1-2) U[T] + 3 m U[T]-2 + m der [U[T] ,T,2] - b; 
In[3] := ATFM[neweq, U, T, 2, c, 1, m, b] 
2 

{a[0] + T a[l] + T a[2], k, c, 1, 0, 0} 
2 2 
4 k 2 c 1 2 2 

{ -2k T , k, c, 1, m, 

3 3 m 2 m 
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2 2 4 4 2 

-16 c - 24 c 1 - 9 1 + 16 k m 
} 

12 m 

which shows the solution 

4:P 2c P / 

u{x, y, z, t) = 2fc^ tanh^ i k(x + ly + mz - ct + d) ] . (64) 

3 3m 2m \ ' 

Here c, d, fc, I and m are arbitrary constants, and b becomes 

-16c2 - 24cP - 91^ + IGk^m^ 



12m 



(65) 



Note that b should vanish for sohton solution in which u — > as j??! — > cxd. In this case 
equation (^5|) is reduced to 

3Z^ + 4c = ±4:k^m. (66) 

Substitution ol the choice 3P + 4c = —Ak'^m into the solution(|6^) gives the famihar 

u{x, y, z, t) = 2fc^sech^ ( k{x + ly + mz - ct + d)] . (67) 



sech^ solution 



5. Discussions 

In this paper, we have obtained the exact N soliton solution of the BS equation 
and the travelling solitary wave solution of equation (p5[). These two solutions seem to 
have essentially the same structure as that of the KdV equation. Indeed their spatial 
dependences are described by a new single variable like pjX + Qjy = p'jx' in equation 



( p^ and x+ly+mz = x' in equation (67). However if we consider multi soliton solution 
and multi soliton collision, the extra dimensions plays essential roles and complex the 
analytical proof of N soliton solutions. V soliton is one of such examples. If we remark 
V soliton collision on some spatial axis we see that two solitons in (1 + 1) dimension 
come together and disappear or that two solitons come to birth from nothing. This 
does not occur in the KdV equation. It is worth noting that this latter process occurs 
in the Broer-Kaup equation which is the (1 + 1) dimensional integrable system written 
in the trilinear form |l^ . 
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Our treatment of extension of integrable system to higher dimensions indicates 
some analogy to that of the d dimensional cylindrical KdV equation. The latter 
system is described by 

ut + ^uux + u^^^ + = O7 (68) 

where d = 1,2 and 3 correspond to the KdV, the cylindrical KdV and the spherical 
KdV equations, respectively. The last term is the curvature term. The Painleve test 
indicates that d = \ and 2 cases are integrable but that d — Z case has the movable 
branch point[^ |l^ This is the same situation as the KdV {d = 1), the BS 

equation {d = 2) and new equations [d — i). However, at last at this stage, it is not 
clear whether these resemblances have any deep implication or not. One of our future 
concerns is to construct an integrable system in (3 + 1) dimensions which is reduced 
to the BS equation and to the KP equation in some particular cases. 
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Figure Captions: 

Figure]^ : Time evolution of the one soliton solution u witli pi — 2, qi — —S. 
Figure H : Time evolution of (j) with pi = 2, qi = —3. 
Figure || : 

(a) An example of the two soliton solution with pi = 0.3, p2 = —0.2, qi = —0.15, 
92 = -0.1. 

(b) Potential diagram corresponding to (a) 
Figure ^ : 

(a) An example of the two soliton solution with pi = P2 = 0.3, qi — —0.15, 
92 = 0.1. 

(b) Potential diagram corresponding to (a) 




EPS File lsspou.ps.ps not found 



Figure 1. 




Figure 2. 





Figure 3. (a) 



(b) 



EPS File 2ssvlures.ps.ps not 
found 




Figure 4. (a) 



(b) 



